function [s,ds] = sigmoid(x)
% modified sigmoid funciton where s = 1/(1+e^{-a(x+b)})
    a = 0.1; b=30;
    s = 1./(1+exp(a.*(x-b)));
    ds = a.*s.*( -1 + s );
end